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ABSTRACT 

We report the discovery of an episode of coherent millisecond X-ray pulsation in the neutron star low-mass 
X-ray binary Aql X-l. The episode lasts for slightly more than 150 seconds, during which the pulse frequency 
is consistent with being constant. No X-ray burst or other evidence of thermonuclear burning activity is seen 
in correspondence with the pulsation, which can thus be identified as occurring in the persistent emission. 
The pulsation frequency is 550.27 Hz, very close (0.5 Hz higher) to the maximum reported frequency from 
burst oscillations in this source. Hence we identify this frequency with the neutron star spin frequency. The 
pulsed fraction is strongly energy dependent, ranging from <1% at 3-5 keV to >10% at 16-30keV. We discuss 
possible physical interpretations and their consequences for our understanding of the lack of pulsation in most 
neutron star low-mass X-ray binaries. If interpreted as accretion-powered pulsation, Aql X-l might play a key 
role in understanding the differences between pulsating and non-pulsating sources. 
Subject headings: pulsars: individual (Aql X-l) — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

Accretion-powered millisecond X-ray pulsars (hereinafter 
AMSPs) had been predicted in the earl y 1980s as the pro- 
genitors of mil lisecond radio pulsars (B acker et all Il982t 
lAlpar et aLlll982l) . The first observational indication that neu- 
tron stars in low-mass X-ray binaries (LMXBs) rotate rapidly 
came in 1996 with the discovery of slightly drifting in fre- 
quency millisecond oscil lations during thermonuclea r X-ray 
bursts (for a review see Strohmayer & Bildsten 2006). How- 
ever it was not until 1998 that th e first AMSP was discovered 
( Wij n ands & van der Klisl lT998). Since then, a total of eight 
AMSPs have been foun d out of the > 150 LMXBs known up 
to date dLiu et al.ll2007l) . 

Since the theoretical prediction of the existence of AM- 
SPs was made, the main issue remained to explain the lack 
of pulsation in the persistent X-ray emission of the major- 
ity of LMXBs. In recent decades many theoretical efforts 
have been made to explain this, the main question remain- 
ing whether the pulsation is hidden from the observer or not 
produced at all. At present, the scenarios most often consid- 
ered are: (a) the magnetic field in non-pulsating LMXBs is 
too weak to allow channeling of the matter onto the magnetic 
poles; (b) the magnetic field has comparable strength inside 
most LMXB neutron stars, but in the large majority of them 
it has been "buried" by a ccretion, resulting in a very low sur- 
face magnetic field (e.g., ICumming et al.ll2001h which again 
is too weak to allow channeling of matter. After the accretion 
stops, the magnetic field eventually emerges from the neu- 
tron star surface and assumes its intrinsic value; (c) pulsations 
are produced in all LMXBs, but in most of them they are 
attenuated by a surrounding scatter ing medium that washes 
out the coherent beamed pulsation (Brainerd & Lambl [19871 
iKvlafis & Klimisll 19871; iTitarchuk et al.ll2003l) : (d) the pulsa- 
tions are attenuated by gravitat ional lensing from the neutron 
star (e.g.. lMeszaros et al.|[l988l) . 

From an observational point of view, since the discovery of 
the first AMSP efforts have been focused on finding differ- 
ences between the sources showing pulsation and those that 

1 Astronomical Institute, "Anton Pannekoek", University of Amster- 
dam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands. E-mail: 
casella@science.uva.nl 



do not, in order to test different theoretical hypotheses. Pos- 
sible observed differences so far are the orbital period (which 
is o n average shorter i n AMSPs than in other LMXBs, see 
e.g., iKaaret et al.ll2006l) and the time-averaged accretion rate 
(which is considered to be o n average smalle r in AMSPs than 
in other LMXBs, see e.g.. lGallowavll2006l) . However, al- 
though it is probable that orbital period and time-averaged 
accretion rate play an important role in the determination of 
AMSP properties, the reason for the lack of pulsation in most 
of LMXBs still has to be found. 

The propert ies of the seventh discovered AMSP (H ETE 
J1900.1-2455. lKaaret et al.ll2006l Idallowav et al.ll2007h gave 
new insights to this issue. This source has a much higher in- 
ferred time-averaged accretion rate than in the other AMSPs. 
The pulsation became undetectable after two months from the 
beginning of the outburst, at strong variance with the other 
AMSPs in which pulsations were always detectable until the 
end of the outbursts. This has been interpreted as evidence for 
burying of t he magnetic field by th e accreted material during 
the outburst (iGallowav et al.ll2007l) . 

Transient hig hly coherent pulsations were als o observed 
in 4U 1636-53 (IStrohmayer & Markwardl 120021) during an 
^800 s interval. However, the pulsations were detected at the 
flux maximum of a superburst, hence they were likely nuclear- 
powered. This is different compared to the seven AMSP, 
where pulsations were interpreted in terms of an hotspot re- 
sulting from magnetic channeling of matter onto the neutron 
star surface. 

We are searching the full Rossi X-ray Timing Explorer 
(RXTE) archive data for coherent pulsations (see also 
I Altamirano et alJfeOOTl) . In this letter we report the discovery 
of an episode of coherent millisecond X-ray pulsarions in the 
neutron star binary Aql X- 1 . Th e LMXB recurren t X-ray tran- 
sient Aql X-l is an atoll source dReig et al.ll2.Q00l) wich shows 
kHz qu asi-periodic oscillati ons, X-ray bursts and burst oscil- 
lations (IZhang etal.lll998bl) . An orbit al period of 18.95 hr 
was obtained by optical measurments (IChevalier & Ilovaiskvl 
I199U IWelsh et al.l2000h. The presumed onset of the p ropeller 
effect (llllarionov & Sunvaevl 19751: IStella et al .11 19861) allowed 
Campanaetal. ( 1998J) to estimate a value for the magnetic 
field of 1 — 3 x 10 8 ( 2 ® kpc ) G, consistent with the magnetic 
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FIG. 1. — Power spectrum of the last 128s data interval of observation 
30188-03-05-00 (see Fig. [2}- The power spectrum was obtained from the full 
bandpass event mode data. The frequency resolution is 1/128 Hz. The pulse 
is clearly visible at ~550 Hz. In the inset we plot the 2-60 keV light curve 
folded at the pulsation period. Two cycles are plotted for clarity. The solid 
line shows the best sinusoidal fit. 

field expected (and later on measured) in AMSPs (see also 
IZhangetal.lll998al) . 

We analyzed RXTE archive data of Aql X- 1 and found that 
the source showed coherent pulsation at a frequency near the 
neutron star spin frequency (as inferred from burst oscilla- 
tions) in its persistent emission for 150 seconds during the 
peak of its 1998 outburst. The discovery of this pulsation 
episode may provide new insights on the issue of why some 
LMXBs pulse and some do not, after many years of debate. 

2. DATA ANALYSIS 

We analyzed the whole PCA public archive of Aql X-l, 
searching for coherent pulsation in its persistent emission. A 
total of 363 observations were analysed. In addition to the two 
Standard Modes always present, the PCA data in different ob- 
servations were obtained in two different high-time resolution 
modes (either Event Mode (125/is) or GoodXenon(~ 1/xs)); 
to all data we applied barycenter correction for Earth and 
satellite motion. For each observation we computed fast 
Fourier transforms from 128s data intervals (after filtering for 
the presence of X-ray bursts) and obtained power spectra with 
a Nyquist frequency of 4096 Hz and a frequency resolution of 
1/128 Hz. We searched all power spectra for significant power 
at frequencies close to the expected pulse frequency (~550 
Hz). Only one detecti on was found with a Leahy power of 
120 dLeahv et al.l 1983b . corresponding to a single trial chance 
probability of 9 x 10 -27 at a frequency of 550.27 Hz (see Fig. 

ID- " 

The single trial significance of the pulsation is 11 cr. Once 
we take into account the number of trials (i.e. the number of 
frequency bins times the number of 128 s power spectra we 
searched) the significance becomes 9 a (3 x 10~ 17 chance 
probability). Combined with its frequency being near that of 
the previously observed burst oscillations, we conclude that 
the pulsation is real and intrinsic to Aql X- 1 and represents the 
spin of the neutron star in this source. A full description of the 
total data analysis is beyond the scope of this letter and will 
be the subject of a forthcoming publication. Here we focus 
on the detailed analysis performed of the single observation 
where pulsation was discovered. 

The pulsation was detected at the end of a 1600-second long 
observation starting on 1998 March 10 at 22:28 UT (Obsld: 
30188-03-05-00). All five proportional counter units were ac- 
tive during the whole observation, with a total 2-60 keV count 
rate varying between 4500 and 5500 counts/sec. During this 
observation high-time resolution data were available in Event 
Mode with 125 /is time resolution and 64 energy channels over 
the 2-60 keV instrument bandpass. The light curve of the ob- 
servation is shown in the top panel of Figure [2] 

The high power was detected in the last 128s interval of the 
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FIG. 2. — Top panel: 2-60 light curve with 16 seconds bin size. Bottom 
panels: Amplitude (middle) and phase (bottom) by epoch folding technique 
every 16 seconds, vs. time. Values of amplitude and phase between and 
~1400 s are consistent with those expected from Poissonian noise. 

observation. By using an epoch folding technique we obtain 
a period estimate of 1.8172746 ms. The folded profile was 
well modeled with a sinusoidal component with 1.9±0.2 % 
amplitude (see inset of Fig. [TJ. No second harmonic could be 
detected, with a 95% upper limit on the amplitude of 0.8 %. 
To study the evolution and the properties of the signal we ex- 
tracted a light curve with maximum time resolution, divided it 
into 16-second segments and determined the pulsation ampli- 
tude and phase in each segment by an epoch folding technique 
(see Fig. |2). The pulsation is detected for ~ 150 seconds at 
the end of the observation; its amplitude gradually increases 
and then decreases again over these 150 seconds. To esti- 
mate the accuracy of our period measurement we performed 
a linear fit to the phases in the interval where pulsations are 
detected. Over this interval the phase is consistent with being 
constant. From a 3 a upper limit of 10~ 3 s _1 on the slope, 
we derive a best period measurement of 1.817275(3) ms. The 
short duration of the pulsation did not allow us to obtain any 
useful upper limit on the frequency drift. From the number of 
cycles we were able to phase connect, we can estimate the co- 
herence of the pulsation as Q > 8 x 10 4 . No correction for the 
binary orbital motion was possible, given the small expected 
drift (< 2 mHz in 150 seconds) and the large uncertainity on 
the orbital phase (~80%, applying the orbital solution from 
IWelsh et alJ2000h . 

In order to study the energy dependence of the pulsation we 
extracted light curves with maximum time resolution in dif- 
ferent energy ranges and determined the pulsation amplitude 
and phase in each of them. In Figure [3] we show the ampli- 
tude vs. energy diagram. The pulsations have a strong energy 
dependence, being not detected at low energies (with an up- 
per limit on its amplitude of ~ 1 % below 5 keV) and strong at 
high energies. 

3. DISCUSSION 

We discovered millisecond X-ray pulsations in the persis- 
tent X-ray emission of the LMXB Aql X-l. The pulsation 
is transient, present for only ^150 seconds and has an aver- 
age fractional amplitude (over the full 2-60 keV instrumental 
bandpass) of ~2%, increasing up to >10% at energies above 
~16 keV. It appears and disappears gradually, on a time scale 
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FIG . 3 . — Amplitude of the pulsation as a function of energy. 1 a errors are 
shown. 

of a few tens of seconds, with a 3 a upper limit on its ampli- 
tude of 0.7% during the ^128 seconds before its appearance. 

This pulsating episode appears to be unique in the whole 
RXTE archive of Aql X-l observations. In one third of the ob- 
servations (corresponding to an exposure of 500 ks) the count 
rate is high enough to allow the detection (3-cr, single trial) 
of a similar pulsating episode as the one we observed. Even 
considering only these data (instead of the total analysed ex- 
posure of 1300 ks) we obtain a recurrence rate smaller than 
3 x 10~ 4 . Such extreme rarety in itself is extremely informa- 
tive. Any physical interpretation must not only in fact explain 
the appearance of the pulsation, but also its extremely low oc- 
currence rate. 

The frequency of the pulsation is 550.27 Hz, ^0.53 Hz 
higher than the reported asymptotic frequen cies for burst os- 
cillations in this source (Zha ng et al.lfT998bh . The light curve 
shows no evidence of an X-ray burst or other obvious feature 
within 1500 seconds before the appearance of the pulsation. 
Together with its high coherence, this leads us to conclude 
that the observed pulsation arises from a hot spot spinning 
with the neutron star. This makes Aql X-l the millisecond X- 
ray pulsar with the longest orbital period and with the highest 
time-average accretion rate known at present. 

The detection of this fast appearing pulsation in a LMXBs 
leads once again to the same question: is the pulsation in Aql 
X-l usually hidden from the observer, by some scattering or 
screening medium that very rarely disappears (or strongly re- 
duces its otpical thickness) for ~150 seconds? Or is the pul- 
sation usually absent, the episode we discovered being the re- 
sult of an occasional and rare asymmetry on the neutron star 
surface? 

In the following we analyse different possible scenarios: 

3.1. Permanent pulsation 

If a hot spot is normally present on the neutron star sur- 
face, the question is what made it observable only for ~150 
seconds. Different authors have argued against or in favor 
of the hypothesis that in non-pulsating sources the pulsation 
is washed out by some optically thick scattering media. The 
main issue is the value of the o ptical depth r needed to atten- 
uate the beamed oscillations. iPsaltis & Chakrabartvl (1 19991) 
discuss the apparent similarities between values of r in pul- 
sating and non-pulsa ting sources, and discard this hypothesis. 
Titarc huk et al.l (120031) report spectral fits giving other values 
of r (for different sources), and concl ude that the data do sup- 
port this hypothesis. More recently, Gogus ^tafl (120071) re- 
port new, independent spectral fits and conclude that in non- 
pulsating LMXBs r is not large enough to cause the pulsa- 
tions to disappear (unless the e lectron temperature is very low, 
see also lTitarchuk et a l. 2007). 

This issue appears to be difficult to resolve, mainly because 



of a substantial degeneracy between optical depth and elec- 
tron temperature in the spectral fits. However, the pulsation 
episode in Aql X-l now allows us to provide new insights 
into this issue. Even though, because of their degeneracy, 
the absolute values of spectral parameters can not be prop- 
erly constrained, one would expect them to change when the 
pulsation appears. In particular, the appearance of the pulsa- 
tion, if caused by a temporary (total or local) decrease of the 
optical depth of the scattering media, is expected to be accom- 
panied by a decrease of the comptonization parameter, hence 
by a softening of the energy spectrum. However, our spectral 
fits show that a two-component model (disk plus boundary 
layer) is sufficient to fit the data, without the need for any 
additional comptonized component. Moreover, neither spec- 
tral fitting nor analysis of hardness light curves show any ev- 
idence of spectral variability associated with the appearance 
of the pulsation. Furthermore, the similarities between the 
spectral properties of Aql X-l and many other LMXBs would 
lead one to expect similar episodes of pulsation in most of 
LMXBs. On the other hand, the very hard spectrum of the ob- 
served pulsation (see Fig. [3) suggests that the beamed radia- 
tion is strongly comptonized (for a discussion on the expected 
energy dependence in the case of the screening scenario see 
iFalanga & Titarchukll2007l) . 

An alternative could be to assume a non-isotropic screen- 
ing medium: an occasional and rare hole in the medium or 
a reflecting cloud could have allowed the radiation from the 
hot spot to avoid the screening medium, thus becoming visi- 
ble. Finally, the hot spot itself could have moved becoming 
temporarily visible, either because of a different absorption 
through the line of sight or because of a change in the rota- 
tional of gravitational lensing geometry. We note that all these 
hypotheses could in principle be consistent with the observed 
gradual increase and decrease of the pulsation amplitude, but 
require very special and fine-tuned geometries. 

3.1.1. Estimate of the magnetic field 

Under the standard dipolar magnetic channeling hypothe- 
sis for creating the hot spot, we can give an estimate of the 
magnetic field. The magnetic field must be strong enough 
to locally disrupt the disk flow. The pulsation episode hap- 
pens when the source is in the soft state, close to the peak 
of the outburst. The 2-10 keV flux is ~ 8.8 x 10~ 9 erg 
s" 1 cm" 2 which correspo nds, assuming a distance of 5 kpc 
dJonker & Nelemansl 120041) a nd a bolometric correction of 
1.4 (Migliari & Fender 2006), to a bolometric luminosity of 
~ 3.7 x 10 37 erg s -1 ~ 0.15L E - Given the relatively 
high luminosity, by imposing that the radius at which this 
happens is larger than the neutron star radius we can de- 
rive a stringent lower limit on the magnetic dipole moment 
(IPsaltis & Chakrabartvl fT999h as \i > 0.47 x 10 26 Gem 3 
which corresponds to a magnetic field of ~ 10 s G at the pole. 

This lower limit is c onsistent with the upper l imit of ~ 
1.8 x 10 9 G obtained bv lDi Salvo & Burderil d2003l) from the 
quiescent lu minosity and with the estimate of 2 ~ 6 x 10 8 G 
obtained by Campa na et alj {1998, assuming a corrected dis- 
tance of 5 kpc instead of 2.5 kpc) from the reported onset of 
the propeller effect. Furthermore, our lower limit rules out 
the hypothesis of residual accretion onto the neutron star to 
explain the quiescent X-ray luminosity in Aql X-l, for which 
an upper limit for the ma gnetic field of ~ 4 x 10 7 G had been 
derived (hypothesis al of lDi Salvo & Burderil (120031) ). 

3.2. Transient pulsation 
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The second possibility is that of a temporary asymmetry on 
the neutron star surface. Let us examine possible causes: 

a) Magnetic channeling 

Accretion rate and magnetic field are the two key physical 
quantities playing a role in channeling of matter. A transient 
episode of channeling of matter could arise due to a temporary 
change in the accretion rate (although the X-ray flux, which 
is thought to be a good indicator of M , does not show any 
variability correlated with the appearance of the pulsation) or 
a change in the accretion geometry. 

An alternative is the occurrence of some variability of the 
magnetic field. In the "burying" scenario, timescales for the 
magnetic field to emerge are set by the Ohmi c diffusion time 
and a re at least of the order of 100-1000 yr (Cu mming et alJ 
1200 ll) . obviously far larger than the few tens of seconds ob- 
served in Aql X- 1 . However a local temporary decrease of the 
Ohmic diffusion timescale could in principle let the magnetic 
field emerge. Possible causes could be a temporary local dis- 
ruption of the screening currents, a starquake, or magnetic re- 
connections. Perhaps, a temporary local magnetic loop might 
form, somewhat similar to those observed in the solar corona, 
where hot coronal plasma "populates" different field lines. 
The magnetic field in this loop might be locally strong enough 
to disturb the flow of matter onto the neutron star surface. The 
resulting hot spot (non necessarily on the magnetic pole) cre- 
ates the observed pulsation. However multipolar magnetic 
fields in neutron stars have not been studied in detail, and 
these phenomena have to our knowledge not been theoreti- 
cally predicted. 

b) Nuclear burning 

The hot spot might be associated to a confined nuclear burn- 
ing event. After its onset, nuclear burning is expected to 
spread over the neutron star surface in a very short timescale. 
Hence one would not expect the resulting pulsation to last 
more than a few seconds. However the burst oscillations 
observed in many sources during the decay of X-ray bursts 
strongly challenge this simplified picture. Moreover, the dis- 
covery of a very coherent pulsation during a superburst in 4U 
1636-53 (Stroh maver & M arkwardt 2002) showed that long- 



lived, very coherent pulsations might arise due to nuclear 
burning. However the lack of any other X-ray signature of 
nulcear burning (such as an X-ray burst) associated with the 
pulsation in Aql X-l strongly argues against this hypothesis. 

c) Surface wave 

An asymmetry on the neutron star surface might result 
from an unstable surface mode, as the ones suggested to 
explain burst oscillations (see e.g., iPiro & Bildstenl l2004t 
iNarayan & Cooper! 120071) . Since the amplitudes of sur- 
face modes have bee n predicted to increase with energy 
dPiro & Bild sten 2006), the hard energy dependence observed 
in Aql X-l, more similar to that of burst oscillations than that 
of accretion-powered pulsations, supports this possibility. 

d) Other possibilities 

The extreme rarity of this pulsation episode leads us to ex- 
amine also more exotic physical explanations. Possible ori- 
gins for a transient hot spot could be a self-luminous magnetic 
event on the neutron star surface, as magnetic reconnections 
or annihilations, or a sudden release of potential energy in the 
deep crust during a starquake. 

Finally, we mention the possibility that some of the dis- 
cussed instabilities c ould arise due to a collisio n with an ex- 
ternal body (see e.g.. lColgate & Petschek 1981). The fluence 
of the observed pulsation is ~ 10 35 ergs, which corresponds 
to the free fall kinetic energy of a mass of ~15 km 3 of iron. 
However, the observed gradual increase in the pulsation am- 
plitude would be difficult to explain within this scenario. The 
probability to have such an impact is presumably extremely 
low. On the other hand, the observed event is very rare. 

Whatever is the physical origin of the observed pulsations, 
their discovery implies the possible occurrence of transient 
pulsations in any of the apparently non-pulsating LMXBs. 
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